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ABSTRACT 

Following the detection of strong TeV 7-ray flares from the BL Lac object 
1ES 1959+650 with the Whipple 10 m Cherenkov telescope on May 16 and 17, 
2002, we performed intensive Target of Opportunity (ToO) radio, optical, X-ray 
and TeV 7-ray observations from May 18, 2002 to August 14, 2002. Observa- 
tions with the X-ray telescope RXTE (Rossi X-ray Timing Explorer) and the 
Whipple and HEGRA (High Energy Gamma Ray Astronomy) 7-ray telescopes 
revealed several strong flares, enabling us to sensitively test the X-ray/7-ray flux 
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correlation properties. Although the X-ray and 7-ray fluxes seemed to be cor- 
related in general, we found an "orphan" 7-ray flare that was not accompanied 
by an X-ray flare. While we detected optical flux variability with the Boltwood 
and Abastumani observatories, the data did not give evidence for a correlation 
between the optical flux variability with the observed X-ray and 7-ray flares. 
Within statistical errors of about 0.03 Jy at 14.5 GHz and 0.05 Jy at 4.8 GHz, 
the radio fluxes measured with the University of Michigan Radio Astrophysical 
Observatory (UMRAO) stayed constant throughout the campaign; the mean val- 
ues agreed well with the values measured on May 7 and June 7, 2002 at 4.9 GHz 
and 15 GHz with the Very Large Array (VLA), and, at 4.8 GHz with archival 
flux measurements. After describing in detail the radio, optical, X-ray and 7-ray 
light curves and Spectral Energy Distributions (SEDs) we present initial model- 
ing of the SED with a simple Synchrotron Self-Compton (SSC) model. With the 
addition of another TeV blazar with good broadband data, we consider the set 
of all TeV blazars to begin to look for a connection of the jet properties to the 
properties of the central accreting black hole thought to drive the jet. Remark- 
ably, the temporal and spectral X-ray and 7-ray emission characteristics of TeV 
blazars are very similar, even though the masses estimates of their central black 
holes differ by up to one order of magnitude. 

Subject headings: galaxies: BL Lacertae objects: individual (1ES 1959+650) - 
galaxies: jets — gamma rays: observations 

1. Introduction 

The EGRET (Energetic Gamma Ray Experiment Telescope) detector on board of the 
Compton Gamma-Ray Observatory discovered 100 MeV-~l GeV 7-ray emission from 66 
blazars, mainly from Flat Spectrum Radio Quasars and Unidentified Flat Spectrum Radio 
Sources (Hartman et al. 1999). Ground-based Cherenkov telescopes discovered TeV 7-ray 
emission from 6 blazars, 4 of which are not EGRET sources. The electromagnetic emission 
of these Active Galactic Nuclei (AGNs) is dominated by a non-thermal continuum with 
a low-energy synchrotron component and a high-energy Inverse Compton component (see 
Coppi (1999), Sikora & Madejski (2001), Krawczynski (2003a) for recent reviews). The TeV 
sources all belong to the class of BL Lac objects, blazars with relatively low luminosity but 
with Spectral Energy Distributions (SEDs) that peak at extremely high energies. 

In the case of TeV blazars, the large detection area of Cherenkov telescopes of sev- 
eral times 10 5 m 2 makes it possible to assess 7-ray flux variations on time scales of min- 
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utes. As the keV X-ray and TeV 7-ray emission from these sources is probably pro- 
duced by electrons of overlapping energy ranges as synchrotron and Inverse Compton emis- 
sion, respectively, observations of rapid flux and spectral variability in both bands com- 
plement each other ideally. The observations can thus be used to constrain and, in prin- 
cipal, even over-constrain models. More specifically, the X-ray and TeV 7-ray observa- 
tions yield a measurement of the jet Doppler factor 5j and the jet magnetic field B at 
the jet base. Observations of TeV blazars can thus reveal key information about the 
astrophysics of mass-accretion onto supermassive black holes and the formation of AGN 
jets. Unfortunately, the interpretation of the TeV 7-ray data is not unambiguous owing 
to the highly uncertain extent of extragalactic absorption of TeV 7-rays in pair-production 
processes with photons of the Cosmic Infrared Background (CIB) and the Cosmic Opti- 
cal Background (CIB). Although X-ray and 7-ray observations of TeV-blazars might ulti- 
mately be used to measure the CIB/COB, a considerable number of sources is needed as 
it is difficult to disentangle source physics and CIB/COB absorption for individual sources 
(Bednarek & Protheroe 1999; Coppi & Aharonian 1999; Krawczynski et al. 2002). 

Owing to its hard X-ray synchrotron emission and low redshift (z = 0.047), the BL Lac 
object 1ES 1959+650 had long been considered a prime-candidate TeV 7-ray source (e.g. 
Stecker, De Jager & Salamon 1996, Costamante & Ghisellini 2002). The "Utah Seven Tele- 
scope Array" collaboration reported the detection of TeV 7-ray emission from the source 
with a total statistical significance of 3.9 a (Nishiyama et al. 1999). The average flux mea- 
sured during the 1998 observations was about that from the Crab Nebula. Motivated by 
the X-ray properties, the Telescope Array detection, and a tentative detection of the source 
by the HEGRA Cherenkov telescopes in 2000 and 2001, we proposed pre-approved pointed 
RXTE target of opportunity observations. These observations were to take place imme- 
diately after a predefined increase in the X-ray or gamma-ray activity was detected with 
the RXTE All Sky Monitor (ASM) or the Whipple 10 m Cherenkov telescope. Following 
the detection of a spectacular TeV 7-ray flare on May 17, 2002 with the Whipple 10 m 
telescope by the VERITAS (Very Energetic Radiation Imaging Telescope Array System) col- 
laboration we invoked the pointed RXTE observations as well as simultaneous observations 
in the radio, optical, and TeV 7-ray bands. The Whipple (Holder et al. 2003) and HEGRA 
(Aharonian et al. 2003b) data showed that the 7-ray flux was strongest during the first 20 
days of observations with peak fluxes of between 4 and 5 Crab units; subsequently, the flare 
amplitude decreased slowly. Following Mrk 421 (z = 0.031) and Mrk 501 (z = 0.034), 
1ES 1959+650 is now the third TeV 7-ray blazar with a high-state flux much stronger than 
that from the Crab Nebula, allowing us to measure the 7-ray lightcurve on a time scale of 
a couple of minutes, and to take energy spectra with good photon statistics on a nightly 
basis. Since the discovery of the first TeV blazar Mrk 421 in 1992 (Punch et al. 1992), the 
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number of well established blazars has now grown to 6 (see Table 1). Fig. 1 shows the 2-12 
keV flux from these 6 sources as measured in the years 1996 to 2003 with the RXTE ASM. 
For Mrk 421, Mrk 501, 1ES 1959+650, and PKS 2155-304 long flaring phases extending over 
several weeks can be recognized. While Mrk 421, 1ES 1959+650, and PKS 2155-304 flare 
frequently, Mrk 501 flared in 1997, but showed only modest fluxes thereafter. The prolonged 
flaring phases offer ideal opportunities to study these objects with high photon statistics. 

In this paper, we discuss the results of the 2002 multiwavelength campaign on 1ES 1959+650. 
We present new radio, optical and RXTE X-ray data taken between May 16, 2002 and Au- 
gust 14, 2002, and combine these data with the already published Whipple and HEGRA 
TeV 7-ray data. In Sect. 2 we present the data sets and the data reduction methods. In 
Sect. 3 we give an overview of the combined light curves, and in Sect. 4 we scrutinize certain 
episodes of the light curves in more detail. After discussing the flux correlations in different 
energy bands in Sect. 5, we present the radio to 7-ray SEDs of 1ES 1959+650 and show 
results of initial modeling with the data in Sect. 6. With the addition of another TeV blazar 
with good broadband data, we consider the set of all TeV blazars to begin to look for a 
connection of the jet properties to the properties of the central engine in Sect. 7. We discuss 
the implications of our observations in Sect. 8. 

We use the following cosmological parameters H = h x 100 km s -1 Mpc -1 with h = 
0.65, Qm = 0.3, and Q\ = 0.7. The redshift of 1ES 1959+650 translates into a luminosity 
distance of 229.5 Mpc. Errors on the best-fit results of x 2_ fits to the RXTE data are given 
on the 90% confidence level. All other errors are quoted on the 1 a confidence level. 

2. Data Sets and Data Reduction 

2.1. Radio Observations 

We used the University of Michigan 26-meter paraboloid to monitor 1ES 1959+650 at 
4.8 GHz and 14.5 GHz between May 5 and August 9, 2002. Each observation consisted of 
a series of ON-OFF measurements taken over a 30-40 minute time period. All observations 
were made within a total hour angle range of about 5 hours centered on the meridian. The 
calibration and reduction procedures have been described in (Aller et al. 1985). Some daily 
observations were averaged to improve the signal-to-noise ratio. 

Additional flux density measurements were made with the Very Large Array (VLA) of 
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the National Radio Astronomy Observatory 15 at frequencies of 43.315 GHz, 22.485 GHz, 
14.965 GHz, 8.435 GHz, and 4.885 GHz on May 7 and June 7, 2002, in snap-shot mode 
(single scans). Observations were made in A- array configuration on May 7 and B- Array 
configuration June 7. Pointing checks were incorporated in an effort to keep the program 
sources near the centers of the primary antenna beams. The data were reduced within the 
A1PS software package supplied by NRAO, following the standard procedures outlined in 
the AIPS Cookbook. Flux density calibration was accomplished through observations of the 
source 3C 286. Some of the flux density measurements at the highest three frequencies were 
discarded owing to erratic variations among antenna pairs. Several secondary calibration 
sources were used for checks on the final flux density scale. 



2.2. Optical Observations 

We present two optical data sets. One was taken with the 0.4 m telescope at Boltwood 
Observatory (Stittsville, Ontario) between May 18 and August 14, 2002, using V, R, and I 
broadband filters. The aperture photometry was performed with custom software and used 
the comparison star 4 from Villata et al. (1998). Data points were obtained from averaging 
over between 4 and 6 2-minute exposures. Relative V and R band magnitudes were converted 
to absolute magnitudes using the published absolute magnitudes from Villata et al. (1998). 
We are not aware of a published measurement of the absolute I magnitude of star 4 and we 
give the results only as relative magnitudes mag(lES 1959+650)-mag(star 4). The typical 
statistical error on the relative photometry of each data point is 0.02 mag. The absolute 
photometry has an additional error of 0.03 mag. 

The other data set was taken with the 0.7 m telescope at the Abastumani Observatory 
in Georgia from May 19 to July 12, 2002, using an R filter for all observations. The frames 
were reduced using DA0PH0T II. The absolute magnitude of 1ES 1959+650 was determined 
by comparison with the standard stars 4, 6, and 7 from Villata et al. (1998). In 20 nights, 192 
measurements of 5 minutes exposure time were taken. The statistical error on the relative 
photometry is 0.1 mag. The absolute photometry has an additional error of 0.05 mag. For 
both data sets we did not attempt to remove the light contribution from the host galaxy. 



15 Thc NRAO is a facility of the National Science Foundation operated under cooperative agreement by 
Associated Universities, Inc. 
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2.3. X-Ray Observations 

The X-ray analysis was based on the 3-25 keV data from the Proportional Counter 
Array (PCA; Jahoda et al. 1996) on board the RXTE satellite. Standard-2 mode PCA data 
gathered with the top layer of the operational Proportional Counter Units (PCUs) were 
analyzed. The number of PCUs operational during a pointing varied between 2 and 4. We 
did not use the 15-250 keV data from the High-Energy X-ray Timing Experiment HEXTE 
(Rothschild et al. 1998) owing to their poor signal to noise ratio. 

After applying the standard screening criteria and removing by hand abnormal data 
spikes, the net exposure in each Good Time Interval ranged from 160 sees to 4.43 ksecs (see 
Table 2). Spectra and lightcurves were extracted with FT00LS v5.1A. Background models 
were generated with the tool pcabackest, based on the RXTE GOF calibration files for a 
"bright" source with more than 40 counts/sec. Comparison of the background models and 
the data at energies above 30 keV showed that the model underestimated the background by 
10%. We corrected for this shortcoming by scaling the background model with a correction 
factor of 1.1. Response matrices for the PCA data were created with the script pcarsp v. 7. 11. 

The spectral analysis was performed with the Sherpa v.2.2.1 package. A galactic neutral 
hydrogen column density of 1.027 x 10 21 cm~ 2 was used for all observations. Since the 
analysis is restricted to the energy region above 3 keV the hydrogen column density has 
only a very minor influence on the estimated model parameters. Single power-law models 
resulted in statistically acceptable fits for all data sets. 

2.4. Gamma-Ray Observations 

1ES 1959+650 was monitored on a regular basis as part of the BL Lac program at the 
Whipple Cerenkov telescope during the 2001/2002 observing season; it was during these 
observations that 1ES 1959+650 was seen to go into an active state. Following the detection 
of strong flares on May 16 and 17, 2002, we coordinated simultaneous observations of 1ES 
1959+650 with the Whipple and HEGRA Cherenkov telescopes. The observations with the 
Whipple 10 m Cherenkov Telescope began on May 16, 2002 and ended on July 8, 2002 
(Holder et al. 2003). The total data set consists of 39.3 hrs of ON-source data, together 
with 7.6 hrs of OFF-source data for background comparison. The Whipple telescope is 
located in southern Arizona, USA, on Mt. Hopkins and is part of the Whipple observatory. 
At this latitude, 1ES 1959+650 culminates at a zenith angle of 33.5°, and so the data were 
necessarily taken at large zenith angles, between 33.5° and 53.5°. The data were corrected 
for large zenith angles and for a temporary reduction of the telescope detection efficiency 
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using the method of LeBohec and Holder (2003) which involves measuring the response of the 
telescope to cosmic rays. While correcting the 7-ray detection rates for the reduced telescope 
sensitivity is straightforward, energy spectra can not be determined with the standard tools 
and further studies of the Whipple energy spectra are underway. The peak energy 16 lies at 
about 600 GeV for the majority of observations. 

Motivated by the HEGRA detection of the source in 2000 and 2001, as well as by the 
strong flaring activity in May, 2002, the HEGRA system of five Cherenkov telescopes (Canary 
Island La Palma) regularly monitored 1ES 1959+650 in 2002. A total of 89.6 hrs of data 
were taken during moonless nights from May 18 to September 11, 2002 (Aharonian et al. 
2003b). Typically, each night comprises about 1 hour of observation time around the object's 
culmination. Owing to the declination of 1ES 1959+650, the object could only be observed 
at zenith angles above 35.9° leading to a mean peak energy of 1.4 TeV. All observations were 
carried out in the so-called wobble mode allowing for a simultaneous measurement of the 
background rate induced by charged cosmic rays. The HEGRA collaboration determined the 
differential 1.3 TeV - 12.6 TeV energy spectrum dN/dE = N (E/l TeV)" r of 1ES 1959+650 
for a high-flux data set and a low-flux data set. The high-flux data set used all 2002 data for 
which the diurnal integral flux above 2 TeV surpassed that from the Crab and gave Nq = 
(7.4+1. 3 stat ±0.9 syst ) 10" 11 photons cm" 2 s" 1 TeV" 1 , and T = (2.83+0. 14 stat ±0.08 syst ). The 
low-flux data set used all 2000-2002 data for which the diurnal integral flux above 2 TeV 
was less than 0.5 Crab units and gave N = (7.8+1. 5 sta t±1.0 syst ) 10" 12 photons cm" 2 s" 1 
TeV" 1 , and T = (3.18±0.17 stat ±0.08 syst ). 

In the following, we quote integral 7-ray flux in Crab units above energy thresholds of 
600 GeV and 2 TeV, for the Whipple and HEGRA data points, respectively. In the case 
of HEGRA, the analysis threshold has been chosen well above the peak energy to minimize 
systematic uncertainties in the region of the trigger threshold. The normalization of the 
fluxes in Crab units renders the results largely independent of Monte Carlo simulations. The 
drawback of the method is that different energy thresholds can introduce flux offsets if the 
source energy spectrum deviates from the Crab energy spectrum. Based on the HEGRA 
results on the correlation of the 7-ray flux level and 7-ray photon index, we estimate that 
these offsets are smaller than 20% for > 600 GeV flux levels on the order of 1 Crab and 
higher, and smaller than a factor of 2 for flux levels well below 1 Crab. Based on the 
Whipple measurement of the energy spectrum from the Crab nebula (Hillas et al. 1998), 
a flux of 1 Crab corresponds to a differential 1 TeV flux of (3.20+0. 17 sta t±0.6 syst ) xlO" 11 
photons cm" 2 s" 1 TeV" 1 and to a vF u flux of (5.12±0.27 stat ±0.96 syst ) xlO" 11 ergs cm" 2 s" 1 . 



16 The peak energy is denned as the energy at which the differential 7-ray detection rate peaks, assuming 
a source with the same 7-ray spectrum as the Crab Nebula. 
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3. Results of the Multiwavelength Campaign 

Figure 2 shows from top to bottom the integral TeV flux, the X-ray flux at 10 keV, the 
3-20 keV X-ray photon index, the V, R, and I band optical data, and the 14.5 GHz and 4.8 
GHz radio data. The TeV 7-ray data (Fig. 2a) show several strong flares during the first 20 
days of the campaign with a flux surpassing 2 Crab units on May 17-20 (MJD 52411-52414) 
and again roughly two weeks later on June 4 (MJD 52429). Subsequently, the flux leveled off 
to about 0.3 Crab units with the exception of two flares on July 11-12 (MJD 52466-52467) 
and July 14-15 (MJD 52469-52470) with a flux between 1 and 1.5 Crab units. Holder et 
al. (2002) studied the fastest 7-ray flux variability time scales based on the Whipple data 
and found a rapid flux increase with an e-folding time of 10 hrs. The large "gaps" in the 
7-ray lightcurves originate from the fact that the Cherenkov telescopes are operated during 
moonless nights only. 

The 10 keV X-ray flux (Fig. 2b, Table 2) was strongest on May 18-20 (MJD 52412- 
52414). It slowly decreased by a factor of 18.7 from the maximum on May 20 to a minimum 
on June 17 (MJD 52442). As we will discuss in more detail in the next section, the TeV 
7-ray and X-ray fluxes seem to be correlated, with the notable exception of an "orphan" 
TeV 7-ray flare on June 4, 2002 (MJD 52429) that is not associated with increased X-ray 
activity. From July 17 (MJD 52469) until the end of the campaign, the X-ray flux stayed at 
a consistently high level: a factor of 1.7 below the maximum flux observed at the beginning 
of the campaign, and a factor of 11.5 above the minimum flux measured on June 17 (MJD 
52442). The TeV emission level during this "X-ray plateau state" is about a factor of 2 lower 
than at similar X-ray flux levels earlier in the campaign. 

We analyzed the X-ray flux variability time scale by computing the e-folding times from 
the flux changes between observations: t — At / A/nF(10keV) with At being the time 
difference between two observations and A/nF(10keV) is the difference of the logarithms 
of the 10 keV fluxes. The shortest e-folding times are given in Table 3. We detected faster 
flux increases than flux decreases: the fastest flux increase has an e-folding time of ~5.9 hrs; 
the fastest flux decrease has an e-folding time of ~15.2 hrs. 

The 3-25 keV photon index T (dN/dE oc E~ r ) (Fig. 2c, Table 2) varies between 1.6 
and 2.4. The X-ray photon index and the X-ray flux are clearly correlated, higher flux 
corresponding to harder energy spectrum. Values well below and well above the value of 
T = 2 show that the low-energy (presumably synchrotron) component sometimes peaked 
above 10 keV, and sometimes below 10 keV. We searched for rapid spectral changes by 
analyzing photon index variations between RXTE observations, see Table 4. The photon 
index V decreased (spectral hardening) by up to 0.09 hr _1 and increased (spectral softening) 
by up to 0.04 fir -1 . As a consequence of synchrotron cooling which is more efficient at higher 
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energies, leptonic models predict that the X-ray emission is harder during the rising phase 
of a flare than during its decaying phase (Kardashev 1962). Careful inspection of the X-ray 
lightcurve and the photon indices does not show evidence for such a behavior. While the 
detection of this effect would impose a constraint on the jet magnetic field and the Doppler 
factor, the non-detection allows large regions in the <5j - B plane (Krawczynski et al. 2002). 

The V,R, and I band optical data (Fig. 2d-f) show flux variations of about 0.1 mag 
on typical time scales of about 10 days. Remarkably, the mean optical brightness increased 
from the first 4 weeks to the last 2 weeks of the campaign by about 0.1 mag in all three 
optical bands. Both, the optical and the X-ray fluxes increased slowly during the campaign. 
Apart from this joint slow flux increase, we did not find any evidence for a correlation 
between the optical and the X-ray or the TeV 7-ray fluxes. We searched for optical intra- 
day flux variability by fitting models to the data of individual days. Although we performed 
very long observations of up to ~7 hrs per night with small statistical errors of 0.02 mag 
per 10 minute exposure time, the reduced chi-square values did not show any evidence for 
statistically significant intra-day flux variability. 

Based on the diurnal brightness averages in the three bands, we computed fastest rise 
and decay e-folding times of 0.07 mag day -1 and 0.03 mag day™ 1 , respectively. Within the 
statistical errors, the V-R and V-I colors stay constant throughout the full campaign. 

The 14.5 GHz and 4.8 GHz radio data taken with the UMRAO (Fig. 2g-h) do not show 
significant flux variations. A fit of a constant flux level to the 14.5 GHz data gives a mean of 
0.174±0.004 Jy with a chi-square value of 24.5 for 21 degrees of freedom (chance probability 
of 27%). The mean flux is consistent with the 14.965 GHz flux of 0.18±0.01 Jy measured 
with the VLA on May 7. A fit of a constant flux level to the 4.8 GHz UMRAO data gives 
a mean of 0.254±0.016 Jy with a chi-square value of 7.1 for 7 degrees of freedom (chance 
probability of 42%). The mean 4.8 GHz flux is compatible with the 4.885 GHz flux measured 
on June 7 with the VLA of 0.23±0.01 Jy, and with 4.85 GHz values of 0.253±0.023 Jy 
and 0.246±0.037 Jy reported by Gregory k Condon (1991) and Becker & White (1991), 
respectively. 

4. Detailed Lightcurves 

In this section we discuss the lightcurves in more detail by dividing the data into 4 
epochs (Epoch 1: MJD 52410-52419, Epoch 2: MJD 52420-52445, Epoch 3: MJD 52460- 
52474, Epoch 4: MJD 52486-52500). In the following figures, we show only the observational 
bands where a substantial number of data points were recorded. 
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Fig. 3 shows the data from Epoch 1 (May 16-25, MJD 52410-52419). Note that the X- 
ray observations started on May 18, UTC 3:26, within less than 24 hrs of the initial detection 
of strong 7-ray flaring activity from 1ES 1959+650 with the Whipple 10 m telescope. The 
7-ray and X-ray fluxes seem to be correlated, both showing a strong flux increase on May 
18 (MJD 52412) and a strong flux decrease on May 21 (MJD 52415). From May 19 to May 
20, the X-ray flux increases by 20% without a similar increase in the 7-ray band. While the 
source brightened in the time interval May 19-25 (MJD 52413.27-52419.40) by about 0.1 
mag in all three optical bands, the X-ray flux decreased by a factor of ~3 over the same time 
interval. 

The data from Epoch 2 (May 26-June 21, MJD 52420-52446) are shown in Fig. 4. 
Except for one strong 7-ray flare, the 7-ray flux stayed well below 2 Crab units. The X-ray 
flux decayed slowly, and the optical brightnesses in the three bands meandered around their 
mean values by 0.05 mag. The most interesting feature of the full observation campaign is 
the strong "orphan" 7-ray flare on June 4 (MJD 52429.308-52429.362) which Fig. 5 shows in 
more detail. While HEGRA measured a low flux of 0.26±0.21 Crab units on MJD 52429.106, 
the Whipple observation 5 hrs later revealed a high flux of 4 Crab units. The X-ray flux 
(measured at the same time as the 7-ray data) did not show any sign of an increased activity: 
the 10 keV flux stayed constant and later even decreased compared to the observation taken 
5 hrs earlier. Similarly, the X-ray photon index and the optical magnitudes do not show any 
irregularity during the 7-ray flare. 

The results from Epoch 3 (July 5-19, MJD 52460-52474) are presented in Fig. 6. The 
7-ray and X-ray fluxes show a very similar development with joint flux minimums on July 13 
(MJD 52468) and July 18 (MJD 52473), and a joint flux maximum on July 15 (MJD 52470). 
As during the full campaign, the X-ray flux and X-ray photon index are tightly correlated. 

The data from Epoch 4 (July 31-August 14, MJD 52486-52500) are shown in Fig. 7. 
While statistical errors hamper the interpretation of the 7-ray data, the X-ray flux varied 
by 50% and the optical flux by 0.1 mag. 

5. Flux Correlations in Different Energy Bands and X-Ray 
Hardness Intensity Correlation 

The correlation between simultaneously measured 7-ray and X-ray fluxes during the 
full campaign is shown in Fig. 8; even though the fluxes seem to be correlated in general, 
the orphan flare clearly deteriorates the quality of the correlation. During the observation 
campaign, the 3-25 keV X-ray photon index and the 10 keV flux were tightly correlated 
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(Fig. 9). Higher flux levels are accompanied by harder energy spectra, as is typical for BL 
Lac objects. The photon-index/flux correlation shows some slow evolution during the multi- 
wavelength campaign with some exceptionally hard energy spectra recorded during Epoch 3, 
around July 15 (MJD 52470). In the hardness-intensity plane we did not detect clockwise 
or anti-clockwise loops during flares. Such loops are expected to occur as a consequence of 
diffusive particle acceleration at strong shocks and synchrotron cooling of the radiating par- 
ticles (electrons or protons) (Takahashi et al. 1996; Kirk & Mastichiadis 1999). The sparse 
observational sampling might be responsible for our non-detection. 



6. Spectral Energy Distribution and SSC Modeling 

The X-ray and 7-ray emission from TeV-Blazars are commonly attributed to the SSC 
mechanism in which a population of high-energy electrons emits synchrotron radiation, fol- 
lowed by Inverse Compton scattering of synchrotron photons to TeV energies. 

In Fig. 10 we show the radio to 7-ray SED of 1ES 1959+650 together with a simple 
one-zone SSC model. The model (Krawczynski 2003b, see Inoue & Takahara (1996), and 
Kataoka et al. (1999) for similar codes) assumes a spherical emission volume of radius R, 
that moves with bulk Lorentz factor T toward the observer. The radiation is Doppler shifted 
by the Doppler factor 

5j = [r(i-/3cos(fl))]-\ (i) 

with (3 the bulk velocity of the plasma in units of the speed of light, and 9 the angle between 
jet axis and the line of sight in the observer frame. The emission volume is filled with an 
isotropic electron population and a randomly oriented magnetic field B. We assume that 
the energy spectrum of the electrons in the jet frame can be described by a broken power 
law with low-energy (E min to E\,) and high-energy (E^, to E max ) indices p± = 2 and P2 = 3, 
respectively (pi from dN/dj oc 7 _Pi , E is the electron energy in the jet frame). Motivated 
by the similar SEDs and flux variability time scales of Mrk 501 in 1997 and 1ES 1959+650 in 
2002, we chose parameter values similar to those inferred for Mrk 501 from time dependent 
modeling of 1997 X-ray and 7-ray data (Krawczynski et al. 2002). The dotted line shows the 
model prediction before taking into account extragalactic extinction and the solid line shows 
the SED modified by intergalactic extinction as predicted by a CIB/COB model with a 
reasonable shape. We choose the CIB/COB model of Kneiske et al. (2002), see e.g. Primack 
et al. (2001) and de Jager & Stecker (2002) for alternative detailed model calculations. The 
parameter values for all subsequent models are given in the respective figure captions. 

While the model shown in Fig. 10 gives a satisfactory fit to the X-ray to 7-ray data, 
it under-predicts the radio and optical fluxes. The model thus suggests that the low-energy 
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radio to optical radiation is dominated by emission from other regions than those that emit 
the bulk of the X-rays and 7-rays. This finding is consistent with the fact that we found much 
less flux variability in the radio and optical bands than in the X-ray and 7-ray bands. The 
Inverse Compton SED corrected for extragalactic absorption peaks in our model at 1.7 TeV. 
Between 100 GeV and 400 GeV the CIB/COB model predicts a characteristic sharp turnover. 
The next generation Cherenkov telescopes CANGAROO III, H.E.S.S., MAGIC, and Whipple 
should be able to measure such sharp turnovers in blazar energy spectra. 

We explored several ways to produce the orphan 7-ray flare in the framework of SSC 
models. Given the observed RXTE energy spectrum and our choice of model parameters, 
it is not possible to produce an orphan 7-ray flare by moving the high-energy cutoff of 
accelerated electrons to higher energies (Fig. 11). The reason for this behavior is that high- 
energy electrons that emit synchrotron radiation above the RXTE energy range, emit Inverse 
Compton 7-rays at energies above those sampled by the observations (above ~10 TeV). The 
additional photons show up at energies above ~10 TeV. Extragalactic extinction reduces the 
flux above 10 TeV already by so much that it is not shown in the figure. 

Adding a low-energy electron population (Fig. 12, left panel) succeeds in producing an 
orphan 7-ray flare and predicts an extremely steep 7-ray energy spectrum. However, the 
model needs careful fine-tuning, as the density of low-energy electrons is constrained by 
the optical measurements. Studies of the TeV 7-ray energy spectrum during the flare are 
underway to test the prediction of a steep spectrum. A more natural way to explain the 
flare is to postulate a second, dense electron population within a small emission region (Fig. 
12, right panel); compared to the region where the quiescent emission comes from, the 1200 
times larger energy density of this electron population and the 5400 times smaller emission 
volume lead to a high Inverse Compton to synchrotron luminosity ratio, and thus to a 7-ray 
flare without a strong X-ray flare. Note that this model does not suffer from a "Compton 
catastrophe" . The optical thickness for internal absorption in 7 TeV + 7seed - ► e+ e ~ pair- 
production processes is well below 1 over the full range of gamma-ray energies covered by 
the TeV observations. 

In Fig. 13 we compare the X-ray and 7-ray energy spectra of 1ES 1959+650 with those 
of the three other TeV blazars with measured TeV energy spectra. The X-ray and 7-ray 
energy spectra of 1ES 1959+650 are very similar to those of Mrk 501. In comparison to 
these two sources, the X-ray spectra of Mrk 421 are softer, while the 7-ray energy spectra 
are similar. The X-ray energy spectrum of H1426+428 seems to be relatively hard. A 
meaningful comparison of the high-energy TeV 7-ray energy spectra of H 1426+428 and the 
other 3 sources is hampered by the highly uncertain extent of extragalactic absorption for 
the high-redshift source H 1426+428. 
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7. Correlation Between Emission Parameters and Black Hole Mass Indicators 

With the addition of another TeV blazar with good broadband data, we consider the 
set of all TeV blazars to begin to look for a connection of the jet properties to the properties 
of the central accreting black hole thought to drive the jet. Ferrarese & Merrit (2000) and 
Gebhardt (2000) discovered a close correlation between the mass of the central black holes, 
M., and the host galaxy's stellar velocity dispersion, a*. The present data on nearby galaxies 
do not show evidence for an intrinsic scatter of the correlation and the upper limit on the 
width of the correlation is 0.4 M.. The correlation is significantly tighter than that of M. 
and the galactic bulge luminosity, L^ g . Based on both correlations Falomo, Kotilainen & 
Treves (2002) and Barth et al. (2003) estimated the black hole masses of several BL Lac 
objects, including 5 of the 6 established TeV-blazars. For 1ES 1959+650, Falomo et al. 
estimated log (M./M ) = 8.12±0.13, using the M.-a* correlation. As can be seen from the 
black hole masses given in Table 1, the black hole of 1ES 1959+650 seems to be the least 
massive of all TeV blazars and is separated by one order of magnitude from the most massive 
one, Mrk 501. 

The black hole mass estimates allow us to explore the correlation between mass and 
parameters describing the jet emission and therefore with the jet properties. In Fig. 14a-f 
we show the correlation between the black hole mass and 6 parameters that characterize 
the jet continuum emission: (i) the luminosity at the peak of the low-energy (synchrotron) 
emission component; (ii) the frequency at which the low-energy SED peaks; (iii) the range of 
observed luminosities at (1+z) TeV; (iv) the 1-5 TeV photon index; (v) the X-ray "flare duty 
cycle", /fl r , defined as the fraction of time during which the RXTE ASM flux exceeds 50% 
of the time averaged flux (see more detailed description below); and, (vi) the range of 7-ray 
e-folding time observed so far. The first two parameters describe the SED of the low-energy 
(synchrotron) component; the second two parameters describe the SED of the high-energy 
(Inverse Compton) component; the last two parameters describe temporal properties of the 
X-ray and the 7-ray emission. 

The 7-ray parameters were corrected for extragalactic extinction based on the CIB/COB 
model of Kneiske et al. (2002). The flare duty cycle was computed from the RXTE ASM 
data taken between 1996 and mid 2003. Binning the data into 28-day bins, we determined 
the fraction of bins where the flux surpassed the mean flux from that source by 50%. We 
determined error bars on these duty cycles with a Monte-Carlo simulation, by re-calculating 
the values 1000 times, modifying each time the flux values according to a Gaussian distri- 
bution with a width given by the experimental flux errors. The choice of the time binning 
changes the values of the flare duty cycle, but does not change the results qualitatively. The 
symbols differentiate the sources; the solid and dashed error bars show the black hole mass 
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estimates based on stellar velocity dispersion measurements from Falomo et al. (2002) and 
Barth et al. (2003), respectively. The dotted error bars show the black hole mass estimates 
from bulge luminosity measurements (also from Falomo et al. 2002). Horizontal error bars 
show the statistical uncertainty on the M.-estimates, and the vertical error bars show the 
ranges of observed values. Differences between parameter ranges can be highly significant 
from a statistical point of view even if the vertical "error bars" span various orders of mag- 
nitudes and exhibit a substantial overlap. The figures do not show clear correlations. The 
only quantity that shows an indication for a correlation is the X-ray flare duty cycle /fl r . 



8. Discussion 

Early SSC modeling of Mrk 421 and Mrk 501 data indicated that simple one-zone SSC 
models were capable of describing a wealth of data satisfactorily (Inoue & Takahara 1996, 
Takahashi et al. 2000, Krawczynski et al. 2001). For Mrk 501 however, detailed time depen- 
dent modeling showed that the very simplest SSC models failed to account for the combined 
broadband X-ray (BeppoSAX, RXTE) and TeV 7-ray data (Krawczynski et al. 2002). In 
order to consistently fit the data from several flares, the authors had to introduce a second 
emission zone, as well as a poorly justified "minimum Lorentz factor of accelerated electrons" 
on the order of 7 m j n = 10 5 and higher. In this paper we presented evidence for an "orphan" 
7-ray flare without X-ray counterpart. Also this finding contradicts the most simple 1-zone 
SSC models. There are several ways to explain the orphan flare: 

• Multiple-Component SSC Models: A high density electron population confined to 
a small emission volume can account for an orphan 7-ray flare (see Sect. 5). "Low duty- 
cycle" fast variability has been observed for a number of sources. A prime example is 
the detection of a strong X-ray flare from Mrk 501 with a doubling time of 6 minutes 
(Catanese & Sambruna 2000). Such observations strongly suggest that indeed small 
regions with high electron densities produce strong and rapid flares. Alternatively, a 
second electron population with a low high-energy cutoff might produce an orphan 
7-ray flare, as mentioned above. However, the corresponding SSC model requires 
fine-tuning of the model parameters. As a third possibility a population of electrons 
with a very hard energy spectrum might produce the gamma-ray flare while it emits 
synchrotron radiation at energies above those sampled with the RXTE. 

• External Compton Models: In External Compton models, the 7-ray flux originates 
from Inverse Compton processes of high-energy electrons with radiation external to 
the jet. Variations of the external photon intensity in the jet frame can cause 7-ray 
flares without lower-energy counterparts. Such variations could have different origins: 
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the external photon flux, e.g. from the accretion disk, could be intrinsically variable. 
Alternatively, the motion of the emission region relative to an external photon reflector 
could result in a time- variable photon flux in the jet frame (Wehrle et al. 1998). In 
External Compton models, the external photon field is highly anisotropic in the jet 
frame, owing to the highly relativistic motion of the jet plasma (r ^> 10). As a 
consequence, the Inverse Compton emission has a narrower beaming angle than the 
synchrotron emission and a slight precession of the jet could cause a large change in 
the TeV flux accompanied by a small change of the X-ray flux. 

• Magnetic Field Aligned along Jet Axis: If the magnetic field in the emission 
region of the orphan flare is aligned with the jet axis and thus with the line of sight, 
the observer would not see the synchrotron flare. The electrons however would scatter 
SSC gamma-rays into our direction and we would thus be able to see the Inverse 
Compton flare. 

• Proton Models: In proton models the low-energy radiation is produced by a popula- 
tion of non-thermal electrons and high-energy radiation by accelerated protons, either 
directly as synchrotron radiation (Aharonian 2000; Miicke et al. 2002), or via a Proton 
Induced Cascade (PIC) (Mannheim 1998). As electron and proton injection rates and 
high-energy cutoffs may vary in a different way with the plasma conditions, proton 
models naturally account for orphan flares. In PIC models, the TeV 7-ray emission 
originates from a thin surface layer of an optically thick pair plasma, while the X- 
ray emission originates from the full emission volume. The model naturally accounts 
for orphan 7-ray flares, as the thin surface layer can produce more rapid flares than 
the larger X-ray emission region. We consider it unlikely that this latter explanation 
applies to the observation of the orphan flare from 1ES 1959+650, as the X-ray and 
7-ray fluxes varied on comparable time scales throughout the rest of the observation 
campaign. 

Our main conclusion from the observation of the orphan 7-ray flare is that it can not be 
explained with conventional one-zone SSC models. 

The black hole mass estimates from stellar velocity dispersion measurements allow us to 
study the connection between the jet emission parameters and the central black hole mass. 
We expect to find correlations as the characteristic length and time scales of the accretion 
system scale with M. (see e.g., Mirabel et al. 1992). Our data however, did not reveal any 
correlations. It is remarkable that 1ES 1959+650 and Mrk 501 show very similar X-ray and 
7-ray energy spectra and flux variation time scales while their black hole masses differ by 
about one order of magnitude. 
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Variations of parameters like jet viewing angle, jet magnetic field, or intensity and en- 
ergy spectrum of the ambient photon field may mask the correlations. Furthermore, our 
correlation plots suffer from the limitations of the observations: flux threshold selection ef- 
fects, limited energy coverage of the observations, and short time over which the data were 
acquired (relative to the lifetime of the jet). Alternatively, the M.-cr* and M,-Lbig corre- 
lations found for nearby galaxies may not hold for blazars, rendering the black hole mass 
estimates used in our analysis inaccurate (Barth et al. 2003). 
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Fig. I. — RXTE All Sky Monitor 2-12 keV lightcurves for the 6 established TeV-blazars. The data 
have been binned to assure a certain minimum signal to noise ratio per point (7 a for Mrk 421 and 
Mrk 501, and 4 a for the other sources). For Mrk 421, Mrk 501, 1ES 1959+650, and PKS 2155-304 
prolonged phases of strong flaring activity can be recognized. 
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Fig. 2. — Results from the 1ES 1959+650 multiwavelength campaign (May 16, 2002 - August 14, 
2002): (a) Whipple (stars) and HEGRA (circles) integral TeV 7-ray fluxes in Crab units above 
600 GeV and 2 TeV, respectively; the Whipple data are binned in 20 min bins, and the HEGRA 
data in diurnal bins; (b) RXTE X-ray flux at 10 keV; (c) RXTE 3-25 keV X-ray photon index; (d) 
absolute V magnitudes (Boltwood); (e) absolute R magnitudes (swiss crosses: Boltwood, circles: 
Abastumani); (f) relative I magnitudes (Boltwood); (g) the 14.5 GHz flux density (UMRAO) and 
(h) the 4.8 GHz flux density (UMRAO). 
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Fig. 3. — 1ES 1959+650 data from Epoch 1 of the campaign (symbols as in Fig. 2). 
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Fig. 4. — 1ES 1959+650 data from Epoch 2 of the campaign (symbols as in Fig. 2). 
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Fig. 5. — 1ES 1959+650 data showing the "orphan" 7-ray flare observed on June 4, 2002 
(symbols as in Fig. 2). 
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Fig. 6. — 1ES 1959+650 data from Epoch 3 of the campaign (symbols as in Fig. 2). 
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Fig. 7. — 1ES 1959+650 data from Epoch 4 of the campaign (symbols as in Fig. 2). 
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Fig. 8. — Correlation between the X-ray flux and the Whipple and HEGRA 7-ray fluxes. 
Epoch 1: full circles, Epoch 2: open circles, Epoch 3: open squares, and Epoch 4: asterisks. 
Only points with a direct overlap of the 7-ray and X-ray observations have been included in 
this graph. 
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Fig. 9. — Correlation between the X-ray flux and the 3-25 keV photon index. Epoch 1: full 
circles, Epoch 2: open circles, Epoch 3: open squares, and Epoch 4: asterisks. 
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Fig. 10.— Radio to 7-ray SEDs of the blazar 1ES 1959+650. The VLA data were taken on May 7 
(filled upward triangles) and June 7, 2002 (filled circles). The optical R-band data (swiss crosses) 
shows the minimum and maximum fluxes detected during the full multiwavelength campaign. Four 
RXTE energy spectra are given (results from the single power-law fits); from top to bottom we show: 
(i) the spectrum during a strong and spectrally hard flare observed on May 20, (ii) an estimate 
of the time averaged spectrum corresponding to the HEGRA "high-state" energy spectrum; (iii) 
the spectrum measured during the "orphan" 7-ray flare on June 4; (iv) the spectrum of the RXTE 
pointing with the steepest photon index (June 14). The open circles show the HEGRA "high-state" 
energy spectrum measured during 6 nights with a >2 TeV integral flux above 1 Crab unit, and 
the diamonds show the HEGRA "low-state" energy spectrum acquired during all 2000-2002 nights 
with an integral flux of less than 0.5 Crab units. An SSC model of the high-state HEGRA data 
and the corresponding high-state RXTE data is shown by the solid line; the dotted line shows the 
model before correction for extragalactic absorption. The model parameters are: 5 S = 20, B = 
0.04 G,R = 5.8xl0 15 cm, log(£ min /eV) = 3.5, log(£ max /eV) = 12.3, log(£ b /eV) = 11.8, Pl = 
2, P2 = 3, electron energy density 0.22 erg cm~ 3 . 
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Fig. 11. — SSC models of the data from the orphan 7-ray flare on June 4, 2002. In addition 
to the data from Fig. 10, the open stars show the TeV flux estimates from 5 hrs before 
and during the orphan 7-ray flare (the X-ray flux stayed at a constant level during the 
flare). The model fits the pre-flare and flare X-ray data, but only the pre-flare 7-ray data. 
The two models computed with different high-energy cutoff of accelerated particles show 
that the additional highest energy electrons mainly produce Inverse Compton emission at 
energies above those sampled by the observations (above ~10 TeV). In the model, the turn 
over of the 7-ray component originates from extragalactic extinction rather than from the 
high-energy cutoff of the electron energy spectrum. The model parameters are: Sj = 20, 
B = 0.04 G, R = 1.4xl0 16 cm, log(£ min /eV) = 3.5, log(£ b /eV) = 11.45, p 1 = 2, p 2 = 
3, electron energy density 0.014 erg cm -3 . Solid line: log(i? max /eV) = 12.2, dotted line: 
log(£' max /eV) = 13.5. All models include the effect of extragalactic absorption. 
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Fig. 12. — Same data as in Fig. 11. In both panels, the solid lines show the SSC model 
that explains the "pre-flare" X-ray and 7-ray emission, and the dotted lines show additional 
emission during the 7-ray flare. All models include the effect of extragalactic absorption. 
In the left panel, the 7-ray flare is produced by an electron population with a rather low 
high-energy cutoff, log(E max /eV) = 11.15 instead of log(-E max /eV) = 12.2. In the right 
panel, a dense electron population confined to a small emission region produces the orphan 
flare. The model parameters for the flare component are as follows, left panel: <5j = 20, B = 
0.04 G, R = 1.4xl0 16 cm, single electron power-law with log(£ , min /eV) = 3.5, log(-E b /eV) = 
log(£' max / eV) = 11.15, pi = 2, electron energy density 0.07 erg cm~ 3 ; right panel: 5j = 20, 
B = 0.04 G, R = 8xl0 14 cm, log(£ min /eV) = 3.5, log(£ max /eV) = 12.2, log(£ b /eV) = 
11.45, p\ = 2, p 2 = 3, electron energy density 17 erg cm~ 3 . The parameters for the quiescent 
emission are the same as in Fig. 11. 
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Fig. 13. — Comparison of the ranges of X-ray and 7-ray energy spectra observed from Mrk 
421 (dotted line), Mrk 501 (dashed line), 1ES 1959+650 (shaded area), and H 1426+428 
(solid line). For Mrk 421, BeppoSAX X-ray data from Fossati et al. (2000) and Whipple and 
HEGRA 7-ray data from Krennrich et al. (2002) and Aharonian et al. (1999c) have been used. 
For Mrk 501, BeppoSAX data from Pian et al. (1998), RXTE data from Krawczynski et al. 
(2000), and HEGRA 7-ray data from Aharonian et al. (1999b, 2001) entered the graphs; the 
different energy coverages of the BeppoSAX (~0. 15-150 keV) and RXTE (3-25 keV) satellites 
resulted in the complex shape of the region of observed X-ray fluxes. For H 1426+428 the 
X-ray data are from (Giommi et al. 2002). and the 7-ray data are from (Petry et al. 2002; 
Djannati-Atai et al. 2002; Aharonian et al. 2003a). The 1ES 1959+650 the X-ray data are 
from this work and the 7-ray data are from (Aharonian et al. 2003b). 
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Fig. 14 . The correlation between the black hole mass estimator and various parameters describing the characteristics of 

the X-ray and 7-ray emission for the five well-established TeV blazars for which black hole masses have been estimated: (a) the 
peak luminosity of the low-energy component; (b) frequency at which the low-energy SED peaks; (c) the range of observed 7-ray 
luminosities at (1+z) TeV; (d) the 1-5 TeV photon index; (e) the flare duty cycle / flr determined from 2-12 keV RXTE ASM 
data, and (f) the range of exponential increase/decay constants observed at 7-ray energies. In all panels, the symbols identify 
the blazars according to the legend given in the first panel; black hole masses from stellar velocity dispersion measurements 
are shown with solid (Barth et al. 2003) and dashed (Falomo et al. 2002) error bars; black hole masses from galactic bulge 
luminosities (Falomo et al.) are shown with dotted error bars. Horizontal error bars are statistical errors in the case of solid 
and dashed lines, for the dotted estimates no statistical errors have been published and we assumed Alog(M, ) = 0.25. The 
vertical error bars show the ranges of observed values. References are given in Table 1; the fastest 7-ray e-folding times for Mrk 
421, Mrk 501, and 1ES 1959+650 are from Gaidos et al. (1996), Aharonian et al. (1999a), and Holder et al. (2002), respectively. 
The available data limit panel (d) to 4 sources and panel (f) to 3 sources. 
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Table 1. Properties of Established TeV-Blazars. 
References: (n) Giommi et al. (2002); (r 2 ) Fossati et al. (2000); (r 3 ) Gaidos et al. (1996); 
(r 4 ) Krennrich et al. (2002); (r 5 ) Aharonian et al. (1999c); (r 6 ) Barth et al. (2003); 
(r 7 ) Falomo et al. (2002); (r 8 ) Pian et al. (1998); (r 9 ) Sambruna et al. (2000); 
(r 10 ) Aharonian et al. (1999a); (r u ) Aharonian et al. (1999b, 2001); (r 12 ) Catanese et al. (1998); 
(n 3 ) this work; (r 14 ) Holder et al. (2002); (r 15 ) Aharonian et al. (2003b); 
(r 16 ) Chadwick et al. (1999); (r 17 ) Horan et al. (2002); (r 18 ) Aharonian et al. (2003a). 
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a Frequency range over which synchrotron peak has been detected. 
b Range of observed synchrotron peak luminosities. 

c Highest (1+z) TcV luminosities observed, including the correction for extragalactic extinction according to Kneiske 
et al. 2002. The number in parenthesis gives the corresponding flux at 1 TeV in (10~ 12 erg cm~ 2 s _1 ). 

d Range of observed 1-5 TeV photon indices, corrected for extragalactic extinction. Numbers in parenthesis give the 
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Mass of central black hole from stellar velocity dispersion. 
f Mass of central black hole from bulge luminosity. 
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Table 2. Results of power-law fits to the 3 keV - 25 keV data 

(Statistical errors only) 
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Table 2 — Continued 
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± 





.017 


1, 


.45 


/ 


45 


0.03 


52472. 


,1680 





.25 


1, 


,012 


± 





.012 


1 


.805 


± 





.017 


0, 


.68 


/ 


45 


0.95 


52473. 


,1562 





.28 


0. 


,963 


± 





.008 


1 


.773 


± 





.012 





.81 


/ 


45 


0.82 


52485. 


,9297 





.54 


1, 


,182 


± 





.005 


1 


.812 


± 





.007 





.99 


/ 


45 


0.49 


52486. 


,9219 





.56 


0, 


,989 


± 





.005 


1 


.875 


± 





.008 





.98 


/ 


45 


0.50 


52487. 


,9688 





.76 


1, 


,089 


± 





.006 


1 


.820 


± 





.007 


0, 


.62 


/ 


45 


0.98 


52488. 


,9570 





.93 


1, 


,094 


± 





,006 


1 


.824 


± 





.008 


0, 


.82 


/ 


45 


0.80 


52489. 


,9453 





.94 


0. 


,925 


± 





.006 


1 


.762 


± 





.009 


1 


.09 


/ 


45 


0.31 


52490. 


.9336 





.94 


1, 


,226 


± 





.007 


1 


.746 


± 





.008 


0, 


.87 


/ 


45 


0.71 
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Table 2 — Continued 



Start MJD t obs a fi 0k cV b T c X ?/d.o.f. d P c c 



52491.9414 


0.46 


1.415 


± 


0.010 


1.759 


± 


0.010 


0.48 


/ 


45 


1.00 


52493.9180 


0.50 


1.003 


± 


0.006 


1.988 


± 


0.008 


0.87 


/ 


45 


0.72 


52494.9727 


0.50 


1.017 


± 


0.008 


1.951 


± 


0.011 


0.72 


/ 


45 


0.92 


52495.9609 


0.53 


1.079 


± 


0.008 


1.973 


± 


0.010 


0.79 


/ 


45 


0.84 


52496.9297 


0.51 


1.362 


± 


0.009 


1.844 


± 


0.010 


0.82 


/ 


45 


0.80 


52497.9297 


0.60 


1.185 


± 


0.008 


1.961 


± 


0.009 


0.76 


/ 


45 


0.88 


52498.9727 


0.56 


1.167 


± 


0.009 


1.981 


± 


0.011 


0.72 


/ 


45 


0.92 



a Duration in hours. 

b 10 keV flux in units of (10~ 3 photons keV^ 1 cm~ 2 s _1 ) 

c 3-25 kcV photon index 

d Reduced x 2 -value and degrees of freedom of the power-law fit 

c Chance probability for larger reduced x 2 -values 
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Table 3. Shortest e-folding times of 10 keV 

flux increases and decreases 



MJDl a MJD2 b At c [hrs] r d [hrs] 



52412. 


.15 


52412. 


36 


4.89 


13.5 ± 0.4 


52412. 


36 


52413. 


.00 


15.40 


-19.9 ± -0.4 


52413. 


.00 


52413. 


.14 


3.49 


5.9 ± 0.2 


52414. 


.47 


52415. 


.14 


16.06 


-19.6 ± -0.3 


52417, 


37 


52417. 


.43 


1.57 


-15.6 ± -1.4 


52429. 


.30 


52429. 


35 


1.10 


-19.6 ± -5.6 


52429, 


35 


52429. 


61 


6.09 


-15.2 ± -0.8 


52432. 


32 


52432. 


.13 


4.49 


-16.5 ± -1.1 


52440. 


.11 


52440. 


.30 


4.68 


10.8 ± 0.7 


52440. 


30 


52440. 


.37 


1.58 


15.3 ± 4.4 


52440. 


37 


52440. 


.43 


1.61 


11.4 ± 2.0 


52444. 


.12 


52444. 


39 


6.26 


-17.2 ± -0.9 



a Centered MJD of first observation 
b Centered MJD of second observation 

c Time difference between observations 

d e-folding time, positive and negative val- 
ues denotes the fastest exponential increase and 
decrease constants, respectively. 
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Table 4. Fastest changes of 3-25 keV photon index. 



MJDl a 


MJD2 b 


At c [hrs] 


Ar/At d [hrs- 1 ] 


52412.36 


52413.00 


15.40 


0.015 ± 0.001 


52413.00 


52413.14 


3.49 


-0.083 ± 0.007 


52431.09 


52431.34 


5.98 


-0.026 ± 0.004 


52440.43 


52441.12 


16.65 


0.015 ± 0.002 


52444.12 


52444.39 


6.26 


0.036 ± 0.004 



a Centered MJD of first observation 
b Centered MJD of second observation 

c Time difference between observations in hours 
d Change in photon index per 1 hr, negative values 
denote spectral hardening 



